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Abstract 

We study the effects of weak point and columnar disorder on the vortex- 
lattice phase transitions in high temperature superconductors. The combined 
effect of thermal fluctuations and of quenched disorder is investigated using 
a simplified cage model. For point disorder we use the mapping to a di- 
rected polymer in a disordered medium in 2+1 dimensions. For columnar 
disorder the problem is mapped into a quantum particle in a harmonic + 
random potential. We use the variational approximation to show that point 
and columnar disorder have opposite effect on the position of the melting line 
as is observed experimentally. For point disorder, replica symmetry breaking 

plays a role at the transition into a vortex glass at low temperatures. 
74.60.Ec,74.60.Ge 
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I. INTRODUCTION 



There is a lot of interest in the physics of high temperature superconductors due to their 
potential technological applications. In particular these materials are of type II and allow 
for partial magnetic flux penetration. Pinning of the magnetic flux lines (FL) by many 
types of disorder is essential to eliminate dissipative losses associated with flux motion. In 
clean materials below the superconducting temperature there exist a 'solid ' phase where 
the vortex lines form a triangular Abrikosov lattice M. In the presence of impurities it was 
suggested H the Abrikosov crystal is replaced by a dislocation -free 'Bragg glass' which 
is also characterized by (quasi-) long range order. This 'solid' can melt due to thermal 
fluctuations or changes in the magnetic field. In particular known observed transitions are 
into a flux liquid at higher temperatures via a first-order melting line (ML) ||, and into a 
vortex glass at low temperature [Q, 0, in the presence of disorder- the so called entanglement 
Une (EL). 

Recently the effect of point and columnar disorder on the position of the melting 
transition has been measured experimentally in the high-T c material Bi2Sr2CaCu20s 
H(BSCCO). Point disorder has been induced by electron irradiation (with 2.5 MeV elec- 
trons), whereas columnar disorder has been induced by heavy ion irradiation (1 GeV Xe or 
0.9 GeV Pb). It turns out that the flux melting transition persists in the presence of either 
type of disorder, but its position shifts depending on the disorder type and strength. 

A significant difference has been observed between the effects of columnar and point 
disorder on the location of the ML. Weak columnar defects stabilize the solid phase with 
respect to the vortex liquid phase and shift the transition to higher fields, whereas point-like 
disorder destabilizes the vortex lattice and shifts the melting transition to lower fields. In this 
paper we attempt to provide a quantitative explanation to this observation. The case of point 
defects has been addressed in a recent paper by Ertas and Nelson using the cage-model 
approach which replaces the effect of vortex-vortex interactions by an harmonic potential 
felt by a single vortex. For columnar disorder the parabolic cage model was introduced 
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by Nelson and Vinokur ||. Here we use an analytic approach to analyze the cage- model 
Hamiltonian vis. the replica method together with the variational approximation. In the 
case of columnar defects our approach relies on our recent analysis of a quantum particle 
in a random potential P,|10| . We compare the effect of the two types of disorder with each 
other and with results of recent experiments. 

Assume that the average magnetic field is aligned along the z-axis which is also the c-axis 
of BSCCO, i.e. perpendicular to the CuO planes. Following EN we describe the Hamilto- 
nian of a single FL whose position is given by a two-component vector r(z) (overhangs are 
neglected) by: 

Here e\ = eo/7 2 is the line tension of the FL, r j 2 = m z /m± is the mass anisotropy, 
eo = ($o/4vrA) 2 , ( $0 is the fluxoid and A is the penetration length), and ji « eo/ a o * s the 
effective spring constant (setting the cage size) due to interactions with neighboring FLs, 
which are at a typical distance of a = \J&o/ B apart. 

For the case of point- disorder, V depends on z and @ 

(V(z,r)V(z',r>)) = Ae 2 e5?\r -r')6(z - z'). (1.2) 

where 

4 2) (r - r') « l/(2< 2 ) exp(-(r - r') 2 /2£ 2 ), (1.3) 

and £ is the vortex core diameter. The dimensionless parameter A is a measure of the 
strength of the disorder. In this case the Hamiltonian is exactly the same as that of a 
directed polymer in a random medium in 2+1 dimensions |IT| . The vortex line represents 
the polymer which is directed along the z-axis (overhangs are neglected). The presence of 
impurities enhances the transverse wandering of the polymer. 

For the case of columnar (or correlated) disorder, V(z, r) = V(r) is independent of z, 
and 
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(V(r)V(r')) = -2/((r - r') 2 /2) = - (1-4) 

In this case we map the problem of a FL in a superconductor into that of a quantum 
particle in a random potential. The partition function of the quantum particle is just like 
the partition sum of the FL, provided one make the identification || 

h—>T, (3fi^L, (1.5) 

Where T is the temperature of the superconductor and L is the system size in the z-direction. 
(3 is the inverse temperature of the quantum particle. We are interested in large fixed L as 
T is varied, which corresponds to high (3 for the quantum particle when % (or alternatively 
the mass of the particle) is varied. The variable z is the so called Trotter time. 
The quantity which measures the transverse excursion of the FL is 

ul{£) ee (\v{z) - r(z + £)f) /2, (1.6) 

The main effect of the harmonic (or cage) potential is to cap the transverse excursions of 
the FL beyond a confinement length £* « O0/7. This length arises by equating the elastic 
energy and the cage potential energy of the FL. Typically after it wanders a distance £* along 
the z-direction the FL is reflected back by the walls of the cage and restarts its transverse 
excursions. The near saturation of Uq(£) at i = £* will become evident from the analytical 
expressions derived in the following sections. We thus define the mean square displacement 
of the flux line by 

u\T)=ul{t). (1.7) 

The location of the melting line is determined by the Lindemann criterion 

u 2 (T m (B)) = c\al (1.8) 

where cl ~ 0.15 — 0.2 is the phenomenological Lindemann constant. This means that when 
the transverse excursion of a section of length rs £*becomes comparable to a finite fraction 
of the interline separation do, the melting of the flux solid occurs. 
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II. THE CASE OF POINT DISORDER 



We start with the case of point disorder that is simpler mathematically. In this case, the 
problem is equivalent to a directed polymer in a combination of a random potential and a 
fixed harmonic potential. The flux line plays the role of the polymer directed along the z-axis. 
The cage potential supplies the harmonic part of the potential and the defects, or impurities, 
the random part. The problem of directed polymers has been investigated extensively in the 
literature in the absence of the harmonic piece. Here we follow the approach of Mezard and 
Parisi (MP) jlT), who used the so called variational (or Hartree) approximation. They set 
up the problem in the presence of a harmonic piece with spring constant /i, but they were 
mainly concerned with the limit of fi — > 0. 

Recall that the Hamiltonian 7i representing the system is given (within the framework of 
the cage model) by eq. ( |1.1| ), together with probability distribution for the random potential 
whose second moment is given by: 

(V(z,r)V(z',r')) = Aeiesf\v~r')6(z-z'). (2.1) 

In order to average over the quenched random potential, the replica method is used. 
After introducing n copies of the fields and averaging over the potential one obtains: 

(Z n ) = [ d[n] . . . d[r n ] exp(-(3H n ) , (2.2) 



with the replicated n-body Hamiltonian given by: 



2 £ 2 , > 

(2.3) 



The variational quadratic Hamiltonian associated with the replica Hamiltonian H n is 
parametrized by: 

l r L 



hn = 75 / dzj^fa ^1 + ^1} 
2 Jo 

1 f L 

~o / dzY\s ab Y a {z) -Y b {z), (2.4) 

2 Jo ~t 



where s ab is the n x n matrix of parameters needed to be determined by the variational 
principles. 

These parameters are fixed by extremizing the variational free-energy given by: 

F = (H n — h n ) hn - i In (/ d[r a ] exp(-^)) . (2.5) 
This free energy is given by: 

F 1 E^f-^M-^/-TrlnGM 



2L 2{3^ uu J-oo2n auy ' 2{3 J 2vr 
where G(u;) is the propagator associated with h n : 



[G = (e, u 2 + fi)5 ah - s ab , (2.7) 

and uj is the 'momentum' variable conjugate to z. Since we are interested in the limit 
L — > oo we can assume that is a continuous variable. Using the formula (for a two 
dimensional vector field r ): 

((r B (*) - r^)) 2 ™)^ = -^[f + <?»M " 2G a6 (^)]J , (2.8) 

we finally arrive at the result: 

— = const. + \ I —(e l cu 2 + /i)VG aa H-i- / — TrlnG(cu) 
2L 2(3 J 2tt 1 2j3J 2tx K ' 

+^Ef P (y^[G aa (u,) + G bb (u)-2G ab (u,)}), (2.9) 
where the function f p is given by: 



Stationarity of the free energy with respect to the parameters s ab gives: 
s ab = 2Pf>A J ^[G aa (co) + G bb (cu) - 2G ab (cu)]), a ^ b 



1 f du, 

~ E *a b . (2.H) 



Here f is the derivative of f p (y) with respect to its argument. 

We consider first the replica symmetric (RS) case, where all the off diagonal elements of 
s a b ar e taken to be equal to each other and their value denoted by s. We denote the value 
of the diagonal elements by Sd- Eq. (|2.11|) implies that in the limit n — > 0, Sd = s. 

In this limit we find: 

GaM = + ^—- 2 , (2.12) 

61 UJ 2 + jl [€i UJ 2 + H) 2 



s = ^f P '(r). (2.13) 



Here we introduced the reduced temperature variable 

r = T/ v ^I. (2.14) 

Using these results we can calculate the mean square displacement ( |1.6| ): 

u 2 (£) = | / ^(1 - cosM))G aa (^) 

= 2T [ ^i^M (1 + _^), (2.15) 

J 2tt ei uj 2 + h K ei u 2 + fi h V ; 

and hence 

u 2 (£)=T(l-e-^)+Ts/ (2») 

x (1 -e- £/e * - (£/£*) e- £/e *), (2.16) 



with 



« 3 A 



T 2ir(£ 2 + r) 

Note that s is positive and independent of u>, and hence the mean square displacement 
ul(£*) is bigger than its value for zero disorder. 

At this point it is convenient to introduce dimensionless variables: 

f = T/(e 0, (2-18) 
B = B£ 2 /%. (2.19) 



In terms of these variables 

H = B e /e, 4 = e/B, (2.20) 

r = £ 2 7 T / y/B. (2.21) 

As parameters for BSCCO we take mean values for those quoted in ref. [|TJ: 

£ = 30A, 

e £ = 1905AT, (corresp. to \ L « 1700A) (2.22) 
7 = 125. 

and also $o — 2.07 x 10 -7 G cm 2 . In the case of the experiments one has 

B < (jf ) 2 , (2.23) 

and hence: 

M 2 (T) / a 2 ~ ( 1 - e" 1 ) + Iy/S(l£) — L- (1 - 2c" 1 ) (2.24) 



2 V V 27r /(7T) ; 



Fig. 1 curve a shows a plot of yuQ(£*)/ao vs. T for A/2n = (curve a) and for A/2tt = 0.2 
(curve b). The value of the magnetic field is taken to be B = 250G. 

Using the Lindemann criterion, eq. ( |1.8|) , we can easily solve for the magnetic field at 
the melting transition: 

B m {T) = ct x (jf ( 1 - e-) + \ (^J ^ (1 - 2e-)) (2.25) 

For T < T cp with 

T cp ^(e e/7)(7A/27r) 1 / 3 (2.26) 
it is necessary to break replica symmetry. This means that the off-diagonal elements of the 



variational matrix s a b are not all equal to each other. MP [11] worked out the equations for 



the replica symmetry breaking (RSB) solution in the limit of \i — > 0, but it is not difficult 
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to extend them to any value of [l. Their final solution is not applicable for the particular 
correlation f p discussed in this paper, hence we will work it out below and in the Appendix. 

When breaking replica symmetry a la Parisi, it is customary to introduce Parisi's pa- 
rameter which is denoted here by u. Thus we put 

Saa = S d , (2.27) 

a j£ b, s ab = s(u), < u < 1. (2.28) 

We have found that a 1-step RSB solution is sufficient in the present case and we parametrize 
it by: 

. , f Sn U < U c , 

s(u) = (2.29) 
Lsi u > u c 

Z = u c ( Sl -s ). (2.30) 

In the following we will used the dimensionless variable s = s£ 2 /eo, and similarly for E. 
In the case of the experiments, the assumption of small \x amounts to the condition (|2.23|) 
which is very well satisfied. For small \x we find (see Appendix for further details): 

u c = f/f cp , (2.31) 

S= ( 7 f cp - 7 f) 2 , (2.32) 
1 B / 7 A\ 

/3 

l T cv = ( ^ ) • (2-34) 



1 B / 7 A\ 

Below T cp the mean square displacement freezes at its value at T cp : 

uliT) ~ u 2 (T cp ), T<T cp (2.35) 

The 'freezing' temperature T cp is about 1.33 larger than the temperature for which the RS 
expression for Uq(T) (see equation ( |2.24| ) has a minimum. The value of the magnetic field 
corresponding to T cp is B m {T cp ) w 2.07($o/^ 2 )(7A/27r) _2 / 3 c| / gives a reasonable agreement 
with the experiments as is evident from Figure 2. 
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In Figure 2 we show a plot of B m (T) vs. T for different values of the disorder. The points 
represent data reported by B. Khaykovitch et al. 0] for various amounts of point disorder 
(we also display points for columnar disorder as will be discussed in the next section). 
The theoretical fit is done using equation ( |2.25| ) for T > T cp and B m (T) = B m (T cp ), for 



T < T cp . Best fit has been obtain by choosing cl = 0.162 and A/2tt = 0.144, 0.208 and 
0.280 respectively. We have chosen the amount of disorder to best fit the temperature 
T cp below which the experimental curves show an apparent change of behavior. This is 
achieved by using eq. ( |2.26|) . The fit associated with the value of A/2tt = 0.144 is the 
one corresponding to the 'as grown' crystal which always has some amount of point defects. 
At high temperatures one observe somewhat larger deviations between the theoretical fits 
and the experimental data. This is due to the proximity to T c « 90K, which affects the 
behavior of the melting line even for pure samples, see discussion in the concluding section. 
The flat part of the curves represent the so called entanglement line, which is believed to be 
a (continuous) transition into the vortex glass. 

III. THE CASE OF COLUMNAR DISORDER 

We consider first the case of columnar disorder. This problem maps into the problem of 
a quantum particle in a random potential, to see this we recall that the density matrix at 
finite temperature (= f3~ l ) of a quantum particle subject to a potential V is given by: 



p(r,r',U) = [dr] exp -- / 

Jr(0)=r ! Tl JO 



™W + EM. + V(r(z)) 



dz\, (3.1) 



with U = f]h.The variable z has dimensions of time and is often referred to as the Trotter 
dimension. This is the same as the partition function of a single flux line, provided one 
makes the identification 

h^T, U = ph^L, m->ei (3.2) 
mentioned in the Introduction. 



10 



In the absence of disorder it is easily obtained from standard quantum mechanics and 
the correspondence (|1.5|), that when L — > oo, 



u 2 (T) = -?=(l- exp(-r sj^i)) = -^=(1 - e" 1 ), (3.3) 



When we turn on disorder we have to solve the problem of a quantum particle in a random 
quenched potential. This problem has been recently solved using the replica method and the 
variational approximation || . Let us review briefly the results of this approach. We apply 
the replica trick in order to carry out the quenched average over the random realizations. 
We consider n-copies of the system, and obtain for the averaged density matrix: 

rr (L)=r n 

p(n--T n ,n--T n ,L) = I " Y[[dr a ]exp{-H n /T}, (3.4) 

Jr a (0)=r a Q=1 

l r L 



I r 



with 



/(y) = -^«p(-^) (3.6) 

In this approximation we chose the best quadratic Hamiltonian parametrized by the 
matrix s a b(z — z')\ 

1 r L 

h n = - dzJ2hr 2 a + vr a ] 

2 Jo 

1 f L f L 

dz dz ' J2 s ab{z ~ z')r a (z) ■ T b (z'). (3.7) 
2T Jo Jo a b 

Here the replica index a = 1 . . . n, and n — > at the end of the calculation. Again, this 
Hamiltonian is determined by stationarity of the variational free energy which is given by 

(F) R IT = (H n - h n ) hn - In J [dv] exp(-h n /T), (3.8) 

The off-diagonal elements of s a b can consistently be taken to be independent of z, whereas 
the diagonal elements are z-dependent. It is more convenient to work in frequency space, 
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where u is the frequency conjugate to z. Uj = (2ir/L)j, with j = 0, ±1, ±2, . . .. Assuming 

replica symmetry, which is valid only for part of the temperature range, we can denote the 

off-diagonal elements of s a b(uj) = (1/T) J L dz e luJZ s a b(z), by s(oj) = s<L,o- Denoting the 

diagonal elements by s d (uj), the variational equations become: 

L ~ , (lT 2T _ 1 \ . , 

s = 2-f c ' — + — £ TTT ^TIS 3 - 9 

~s d (u) = ~s-% f L d( (l-e^)x 



l '\L^ e^ + l-M-'))- (3 ' 10) 
here f c '(y) denotes the derivative of the "dressed" function f c (y) which is obtained in the 
variational scheme from the random potential's correlation function f(y), and in our case is 
given by: 

Uv) = - 9J ^7^- (3-11) 

The full equations, taking into account the possibility of replica-symmetry breaking are 
given in ref. ||. In terms of the variational parameters the function Uq(£) is given by 

*w=gs Jr™ ( ^iU - (3 - 12) 

^ w ,^ tiUJ + fJ, — S d {UJ ) 

This quantity has not been calculated in ref. |f§. There we calculated (r 2 (0)) which does 
not measure correlations along the ^-direction. 

In the limit L —>■ oo we were able to solve the equations analytically to leading order in 



g. In that limit eq. ( |3.10|) becomes (for lu ^ 0) : 



4 ^ 2 r°° 

s d {u) = -fc "(b ) - - / efc(l - cos(^)) 
fi 1 Jo 

x(AW-iM ( 3 - 13 ) 

with 

CM-zrrp. 1 - cos(w) . , , (3.14) 

K = 2tT^ 1 (3.15) 

J-oo 2ir eiuo z + \i — s d {uj) 

12 



One can solve equation ( |3.13|) numerically, but there it is hard to obtain good accuracy at 
high frequencies when the cosine term oscillates strongly. We can get a better approximation 
analytically. We parametrize an approximate solution by: 

s d (co) = Soo + Afi/(eiUJ 2 + a 2 fi), (u ^ 0) (3.16) 

and require that it will obey the correct behavior at low and high frequencies to leading 
order in the strength of the disorder. There are three parameters which are determined by 
Sdioj = 0), s' d {u = 0), and Sd(uj = oo). To leading order in the strength of the disorder we 
can substitute in eq. (|3.13|) , 

bo = r, C (C) = r{l - exp(-|CI ), (3.17) 
with t = T / ^Jti fi. We then find after some algebra 

s - = -I "(r)(l + 1 j[ °dC e'H-^—, + n 1 _ 02 )) 
fi A Jo 1 — ae ^ (1 — ae 

= i/ c ff (r)(4 + / 1 (a)), (3.18) 

where we defined 

a = r/(e 2 + r), (3.19) 

h{a) = 1/(1 - a) - (1/a) log(l - a). (3.20) 

Similarly for small u we find: 



Sd (o;) = -f c "(r)(l + c V / rfC C 2 e^( i - , + 



with 



/i 8 fi Jo ' 1-ae ^ (1- ae~^) 2 

l/ c "(r)(l + ^^-/ 2 («) + ...), (3-21) 



1 00 J- _|_ 1 

/ 2 («) = -E^« fc - (3-22) 

" k=l ft 



From this equation we find for the other two parameters in eq. (|3.16|) 
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a' = M<*)/Ma), (3.23) 

A= _f^V)fM. (3.24) 
(i / 2 (a) 

Notice that Sd(u;) is negative for all u > 0. It interpolates from the value — 4| / c "(t)| ///at 
~ to the value -(4 + /i(a))| f c "(r)| //i at w = oo. Substituting ( ^TB|) in eq. ( gig ) 
and expanding the denominator to leading order in the strength of the disorder, we get : 

u 2 (l) = C (l) = r(l - A/(a 2 - I) 2 / ft) 

x(l - e~ £/r ) + TA/(a(a 2 - l) 2 fi) x 
(1 _ e -^*) + r /( 2/i ) x ( Soo + A/(a 2 - 1)) 
x (l-e-W -(£/t) e~ e/i *). (3.25) 

Recall that this result is valid to first order in the strength of the columnar disorder. 
This expression simplifies significantly under the assumption 

B < ( 7 T ) 2 , (3.26) 

valid in the experiments. In that case 

a ~ 1 - / 7T + • • • , (3.27) 
a 2 ~ 0.3517T / \/S - 0.475 log(Vs" / 7?) + 0.326 + • • ■ , (3.28) 
A/fl ~ o, 35 i|_ _J_ (i _ ( i.07 + 2.351 log(V5~ / tT))V^~ / 7^) + ■ • • , (3.29) 

WA* = -^ Bfrf)^ 1 + {2 ~ l0g(V/ ^ 7 7? ))V ^" / lf) + ' " ' (3 - 30) 
from which we find: 

u 2 (T) / a 2 , ~ 0.632 /^7T - 1.952 f-^j — L- + 4.804 f— ^ —2 — - + •- • . (3.31) 
v ;/ v ' \2nJ ( 7 T) 2 V 27r /(7^) 5/2 

From this equation we can derive the most important result for the location of the melting 

transition: 

B m (f) = 2 -^(l + 3.904c- (JL) JU(1 - ™%*)) . (3.32) 
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In Fig. 1 shows a plot of yw 2 (T)/ao vs. T for g/2ir = 0.025 (curve d). We have chosen 
B = 250G. We see that the disorder tends to align the flux lines along the columnar defects, 
hence decreasing u 2 (T). Technically this happens since Sd{w) is negative. 

For the case of combined columnar and point disorder, it is tempting to combine equa- 
tions ( 2.25 ) and ( |3.32 ) into a single equation for T > T cp : 



9 .n dr 4 l + 3.904c^ 2 f Ir-Ml- 2 ^) 
BUT) = ^ ; UJ(7 ! )2 * • (3.33) 



0.209 f^T^ 



( 7 T)3_ 

Even the naturally grown crystals has some amount of point disorder as discussed above. 

In Fig. 2 we show the modified melting line B m (T) in the presence of columnar disorder 
and a small amount of point disorder ( A/2ti = 0.144, corresponding to the 'as grown' curve 
fit of the last section), as given by eq. ( |3.33| ) with cl = 0.162. We see that the melting line 
shifts towards higher magnetic fields with increasing amounts of columnar disorder. The best 
fit to the experimental results is obtained for g/2ir = 0.01 and g/2n = 0.025 respectively. 

For T <T c k, (e £/7)[45( 2 <I> / (tt 2 ^ 2 B)} 1/6 , there is a solution with RSB. This temperature 
is below the bottom of the range plotted in the figures for columnar disorder. It is thus not 
necessary to include the RSB solution in the plot. 



IV. CONCLUSIONS 

The analytical expressions given in eqs. ( |2.25|) , ( |3.32| ), though quite simple, seem to 
capture the essential feature required to reproduce the position of the melting line. The 
qualitative agreement with experimental results is remarkable, especially the opposite ef- 
fects of point and columnar disorder on the position of the melting line. The 'as grown' 
experimental results are, as expected, corresponding to small amount of point disorder. 

The effect of point disorder is to increase the transverse excursions of the FL which 
seeking the best free-energetic configuration. This increase in the mean square fluctuations 
lowers the melting temperature for a given magnetic field. At low temperature, the entan- 
glement transition is associated in our formalism with RSB, and is a sort of a spin-glass 
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transition in the sense that many minima of the random potential and hence free energy, 
compete with each other. This means that there are many possible deformations of the FL 
which are very close in energy, or free energy. To our knowledge this is the first time the 
transition into the vortex glass is represented as a RSB transition. 

For point disorder, in the limit of infinite cage ( /x — > 0), the variational approximation 



gives a wandering exponent of 1/2 for a random potential with short ranged correlations [11 



whereas simulations give a value of 5/8 ||12|| . This discrepancy does not seem of importance 



with respect to the conclusions obtained in this paper since we always consider the case of 
finite n (which amount to a non-zero magnetic field) and also consider distances of order £* 
in the z-direction. 

For columnar disorder, the z-independence of the random potential tends to reduce the 
transverse fluctuations of the FL, thus shifting the melting transition to higher temperatures 
and magnetic fields. Related to this point is the fact that columnar disorder is much more 
effective in shifting the position of the melting line as compared with point disorder for 
the range of parameters considered here. We have used a much weaker value of correlated 
disorder to achieve a similar or even larger shift of the melting line than for the case of 
point disorder. This is again related to the z-independence of the random potential in the 
columnar case, which help to enhance its effect on the excursions of the vortex lines. 

The experiments show that in the case of columnar disorder the transition into the vortex 
glass seems to be absent. In our model it corresponds to the fact that RSB does not occur 
in the temperature range of relevance to the experiments, but rather at significantly lower 
temperatures for the amount of disorder present. 

Concerning the apparent deviations of the theoretical curves from the experimental points 
at high temperatures, we should point out that our expressions are only valid far from Tc 
which is ~90K for BSCCO. Close to T c the melting line cannot behave as B m oc 1/T 2 since it 
must terminate at B=0, T = T c . This comes about because of the temperature dependence 
of the fundamental constants like the penetration and coherence lengths. Houghton et al. 
13| used a detailed melting theory to obtain a behavior of B m rs Bq(1 — T/T c ) 2 near Tc. 
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Other authors JT^] predict a behavior like B m pa Bq(T c /T — 1), which seems to better fit the 



experimental data for the system under consideration according to ref. ||. Adding such a 
correction factor to our expression for B m (T) can improve the fit at high temperatures. 

Another point we should mention is that our model incorporates the disorder as depend- 
ing on a single parameter like A or g , whereas experimentally there are two parameters 
which are the density of defects, which is varied experimentally, and their individual strength 
over which there is not much control. Our model is valid in the limit of very weak impurities 
which are densely and uniformly distributed in space. This is quite reasonable for the case 
of point disorder, but for columnar disorder the experiments involve a rather low density 
of columnar defects which is smaller than the density of vortices. Thus from the point of 
view of a single vortex the disorder is not uniformly distributed which may account for the 
apparent difference in curvature between theory and experiment. 

We have shown that the cage model together with the variational approximation repro- 
duce the main feature of the experiments. Effects of many body interaction between vortex 
lines which are not taken into account by the effective cage model seem to be of secondary 
importance. Inclusion of such collective effects within the variational formalism remains a 
task for the future. These effects may be responsible for the apparent difference in curvature 
between the experimental and theoretical curves, for the case of columnar disorder. 

We thank David Nelson and Eli Zeldov for discussions. We thank the Weizmann Institute 
for a Michael Visiting Professorship, during which this research has been initiated. 

V. APPENDIX 

In this appendix we give more details of the RSB solutions for point disorder and for 
columnar disorder. We use some of the results found in Appendix II of MP ( ]TT]). Using the 
parametrization given in eqs.( |2.29"D and ( |2.30| ), we find for the free energy: 



F rl-ii„_, / u t 1 — u.. /— — 

— = canst. + ~ 9 VW + £ 
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Sir £ 2 + — — iisr./A o7T £ 2 i r -H— 



where we defined A = Acq£ 3 . Stationarity with respect to £ and w c yields equations for 
these two quantities. Introducing the dimensionless quantity 

£ = ££ 2 /e , (5.2) 

and taking the limit of small /i, these equations become: 

~ u c 7A 1 



7 T 27r (I + 7T /\/s) 2 ' 



(5.3) 
(5.4) 

(7^) 2 ^(1 + 7 T/^) 
Solving these equations we find the solutions given in equations fl2.31|) and (|2.32|) . The 
quantities s and si are given by the equations: 

s = 2/3— r j^^, (5-5) 

si = 2/3— , (5.6) 

from which equation ( p. 33 ) follows in the limit of small /i. 
The mean square displacement Uq(£) is given by 

u 2 (£) = 2Tj ^(1 - co8M)G M M 



-(1 - exp(-£v^)) - L_%- /_^_(i _ exp(-£ V / ( yU + S)/ e/ )) 
+s ^(l - exp(-£ v /^) - (fy^) exp(-£ v /^)). (5.7) 



From this expression we obtain eq.( [2.35 ) for small /i. 



Let us discuss briefly the case of columnar disorder. In this case we showed in ref. ||that 
replica symmetry is broken in a region of the T — h 2 jm phase diagram. In the present paper 
we use the mapping given by equation ( |3.2j ), and in addition we take the limit L — > 00 . 
For the case of short ranged correlation of the potential a 1-step RSB has been found. In 







u < u c 


s(z) 


= 








u > u c 


E 


= u c (si 


- so), 


u c 


= yjL, 





the limit L — > oo we find that the breaking point u c of the one step solution scales as L^ 1 . 
Nevertheless there is a finite contribution to observable quantities as will be discussed below. 
Putting 

(5.8) 
(5.9) 
(5.10) 

we find the following equation for the limit of large L: 

9 ~ 9TY 
s /L = -f c ' (b + ) (5.11) 

T y c n(/i + E) 

si/^ = |/ c / (M (5.12) 

a d (u;) = -E + — =/ c " (6b) - - / *(1 - cos(c^)) 
/x + E i Jo 

x(/ c '(C7 (?))-/c'(6b)), (5-13) 
and E and y c satisfying the equations 

S = ^(/c ' (6b) " /c ' (6b + ), (5-14) 

T y c /i(/i + E) 

= — T 2 E + T 2 (/i + E) log(l + E//i) 

- - 2TE 

// y c /x(/x + E) 

These equations possess a RSB solution (they always posses the solution E = 0) below a 
temperature T c , which for small \x is given by 

1/3 



= b- 1 ". (5.16) 



Just below T c we find to leading order in T c — T : 



E~f*F~^), (5.17) 



y c £_^(i_3±£_jL). (5.18) 
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Figure Captions: 

Figl: Transverse fluctuations vs. temperature in the cage model for fixed B = 250G. 
(a) no disorder (b)point disorder (A/(27r) = 0.2) (c)RSB for point disorder, T < T cp , 
(d)columnar disorder (g/(2n) = 0.025 ) 

Fig. 2: Melting line for different amount of point and columnar disorder. The experi- 
mental points from ref. || are denoted by symbols and the theoretical prediction by lines, 
(a) squares: as grown sample, continuous curve A/(27r) = 0.144. (b) circles: point disorder 
induced by a dose of 3 x 10 18 e~/cm 2 , dashed line A/(27r) = 0.208. (c) triangles: point disor- 
der induced by a dose of 6 x 10 18 e~ /cm 2 , dotted line A/(27r) = 0.280. (d) inverted triangles: 
columnar disorder equivalent to = 50G, dashed-double-dotted line g/(27r) = 0.01 and 
A/(27r) = 0.144. (e) diamonds: columnar disorder equivalent to B^ = 100G, dashed-dotted 
line g/(2n) = 0.025 and A/(2tt) = 0.144. 
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